Advanced optimal extraction for the Spitzer/IRS 



V. Lebouteiller, J. Bernard-Salas, G.C Sloan, and D.J. Barry 

Center for Radiophysics and Space Research, Cornell University, Space Sciences Building, Ithaca, NY 

14853-6801, USA 
vianneyOisc . astro . Cornell . edu 



ABSTRACT 

We present new advances in the spectral extraction of point-like sources adapted to the 
Infrared Spectrograph onboard the Spitzer Space Telescope. For the first time, we created a 
super-sampled point spread function of the low-resolution modules. We describe how to use the 
point spread function to perform optimal extraction of a single source and of multiple sources 
within the slit. We also examine the case of the optimal extraction of one or several sources with 
a complex background. The new algorithms are gathered in a plugin called AdOpt which is part 
of the SMART data analysis software. 

Subject headings: methods: data analysis, techniques: spectroscopic, infrared: general 



1. Introduction 

The ideal spectral extraction algorithm for 
point-like sources yields the maximum signal-to- 
noise ratio (S/N) while at the same time preserv- 
ing the spectrophotometric fidelity. The standard 
extraction method is based on co-adding the flux 
in the cross-dispersion direction within a window 
large enough to contain (most of) the source flux. 
This method is known as "tapered-column" in the 
"Spectroscopy Modeling Analysis and Reduction 
Tool" (SMART0, Higdon et al. 2004) and is equiv- 
alent to the "regular extraction" of the Spitzer 
IRS Custom Extractor (SPICE), provided by the 
Spitzer Science Center (SSC). Although this gen- 
erally produces satisfactory results, the inclusion 
of noisy pixels which do not contain a significant 
fraction of the flux inevitably tends to degrade the 
quality of the extracted spectrum. This is because 
every pixel in the extraction window is given the 
same weight. Moreover, the extraction window is 
part of a pseudo-rectangle which is defined as a 
zone in the detector array where the wavelength is 
uniform (Figure[T|). When a quadrilateral bound- 
ary crosses a pixel, the signal is assumed to be 
evenly distributed within that pixel, which can 



lead to artificial flux variations in the extracted 
spectrum. Because of the interplay of the angled 
spectral trace and the widening extraction aper- 
ture, this error oscillates with wavelength with an 
amplitude which decreases toward longer wave- 
lengths. 

Knowledge of the point spread function (PSF) 
of the instrument can solve these problems, as 
it enables the so-called optimal extraction tech- 
nique which weights the extracted data by the S/N 
of each pixel (Horne 1986). Optimal extraction 
therefore significantly reduces the statistical noise 
in the final spectrum compared to more typical 
extraction algorithms, which weight all the data 
within the window equally. So far, efforts on op- 
timal extraction for the IRS have made use of 
template PSFs or analytical PSFs to fit the cross- 
dispersion profile of the data: 

• Virtually any spectrophotometric standard 
star can be used to derive a template PSF, 
which then allows an empirical estimate of 
the PSF at the reference positions (referred 
to as "nod" positioned). This method is cur- 
rently used by the software SPICE (see Nar- 



1 SMART is available at http://isc.astro.cornell.edu/smart/. 



2 See the Spitzer/IRS observer's manual 
http:/ / ssc. spitzer. caltech. edu/ documents /SOM/ 
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ron et al. 2007). The resulting extraction 
undeniably provides better S/N as compared 
to a tapered column extraction, although the 
corresponding algorithm is sensitive to cross- 
dispersion offsets between a source position 
and the nod position. A simple shift of the 
PSF is not sufficient to acquire the best S/N 
possible, and in some cases, low- frequency 
oscillations can appear in the spectrum due 
to the misalignment. The reasons are inher- 
ent to the data used to compute the PSF 
template since the latter is created in a spe- 
cific observational mode (default positions 
along the slit, default data sampling). 

• Analytical PSFs allow estimating the instru- 
mental profile for any position along the 
slits. The independent effort from the "Core 
to Disks" legacy program (Evans et al. 2003) 
uses such PSFs along with an extended emis- 
sion background which is determined on-the- 
fly (see Lahuis et al. 2007). Analytical PSFs 
can bear some uncertainties due to the lack 
of knowledge of the exact instrumental pro- 
file. 

We have developed a new optimal extraction 
algorithm to be used with either of the low- 
resolution modules of the Infrared Spectrograph 
(IRS, Houck et al. 2004) onboard the Spitzer Space 
Telescope (Werner et al. 2004). The new extrac- 



Fig. 1. — The standard extraction aperture for a 
single wavelength, superimposed on the grid de- 
fined by the pixels on the detector array. The 
quadrilateral represents the extraction window, 
which is a "tapered-column" whose width in- 
creases proportionally with wavelength to account 
for the varying point spread function. The tilt 
of the quadrilateral reflects the fact that the row 
axis of the detector array is not parallel to a line 
of constant wavelength. 



tion algorithm is available via the plugin AdOptjj, 
part of the new release of SMART. In a nutshell, 
an empirical super-sampled PSF has been con- 
structed for each row of the detector array and a 
multi-linear regression algorithm is used to weight 
the pixels and derive the flux from the source. 
The optimal extraction is based on detector rows 
rather than pseudo-rectangles to treat each pixel 
as indivisible and thus avoid uncertainties due 
to the lack of knowledge on the pixel response 
function. The algorithm presented in this pa- 
per produces a considerably higher S/N, up to a 
factor of ~ 2 for faint sources, than the current 
extraction method available in SMART for point- 
like sources ("tapered column"). An additional 
advantage of using a super-sampled PSF is that it 
remains valid anywhere along the aperture in the 
cross-dispersion direction. We found that the op- 
timal extraction method is extremely sensitive to 
offsets between a source position and the nominal 
position. For this reason a new algorithm to locate 
the source in the slit has been implemented, with 
a precision of better than a twentieth of a pixel. It 
is also now possible to extract spectra of spatially 
blended sources, which is crucial when dealing 
with crowded regions, such as stellar clusters or 
nearby galaxies. Such as what can be achieved 
using iterative techniques (Lucy & Walsh 2003), 
the cross-dispersion profile of the data is decom- 
posed into its components, including the spatial 
profiles of any number of sources in the slit and 
the extended background emission. Finally it is 
also possible to extract sources with significant 
offset in the dispersion direction by calculating a 
modified PSF on-the-fiy. 

In the following section the basic steps to con- 
struct the PSF are given. Section[3] details the 
mathematical description of the method and its 
application for extracting multiple sources are de- 
tailed. Section|4] briefly explains some specific ap- 
plications. 

2. Super-sampled PSF 

We have constructed super-sampled PSFs for 
the two low-resolution modules (Short-Low, SL 
and Long-Low, LL, see Tabled]) of the IRS on 
Spitzer. The optimal extraction for the high- 



3 The documentation is available at 

http://isc.astro.corneU.edu/SrnartDoc/SrnartOptimal. 
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Table 1: Main properties of the low-resolution modules of Spitzer/IRS. 



Module 


Order 


A ( fj,m) 


Aperture size (") 


Pixel size (") 


SL 


1 


7.4 - 14.5 


3.7x57 


1.8 


SL 


2 


5.2 - 7.7 


3.6x57 


1.8 


SL 


3 


7.3 - 8.7 


3.6x57 


1.8 


LL 


1 


19.5 - 38.0 


10.7x168 


5.1 


LL 


2 


14.0 - 21.3 


10.5x168 


5.1 


LL 


3 


19.4- 21.7 


10.5x168 


5.1 



resolution modules will be considered in the future 
since these modules do not allow the full sampling 
of the PSF profile. 

In order to build the super-sampled PSF, 
we considered observations of calibration stars 
scanned around reference positions in the cross- 
dispersion direction. We used the basic-calibrated 
data (BCD) product. After cleaning the bad 
pixels using IRSCLEAN^, individual exposures 
(DCEs) were combined, after verifying that the 
pointing remained stable. The source in each im- 
age was then found via the source finder algorithm 
(Sect. [3?6)) which makes use of previous iterations 
of the PSF profile. Finally, the positions of the 
source along the slit were used as inputs in an iter- 
ative reconstruction of the high-resolution spatial 
profile from the under-sampled data. We used 
a simplified version of the regularized image re- 
construction explained in detail by Pinheiro da 
Silva et al. (2006). In short, considering a series of 
/c-undersampled spectra D k and the desired high- 
resolution profile P, the observed spectra take the 
form: 

D k = UW k P + r )k , (1) 

where U is the downsampling matrix, W is the ge- 
ometric transformation matrix, i.e., in our case the 
relative shift with respect to a reference position, 
and rj represents the noise. W is easily calculated 
from the source positions, while U is the matrix 
corresponding to a simple re-scale with no interpo- 
lation. We used the same algorithm described by 
Pinheiro da Silva et al. (2006), with the following 
iterative step: 

P/ +1 = Pj+^iWZuTD^-iwZuTUW^Pj-XiPj 

k 

(2) 



4 IRSCLEAN can be found at http://ssc.spitzer.caltech.edu 



where fi controls the convergence speed, A con- 
trols the regularization, and Pq is a reference spec- 
trum, which can be either the theoretical PSF, or 
an initial guess based on a simple co-addition of 
the aligned data images. 

In practice, we calculated the SL PSF on a ref- 
erence grid with 10 sub-pixels per pixel per actual 
pixel, while the LL grid contains five sub-pixels. 
The number of sub-pixels was chosen to maximize 
the signal in each sub-pixel while allowing the best 
resolution possible on the PSF profile. The shape 
of the PSF is the same within the uncertainties for 
the two nod positions. Table [2] gives the coordi- 
nates for the nod positions and the equations defin- 
ing the spectral trace in each module. The spectral 
trace quantifies how the position of the centroid of 
the PSF shifts with wavelength (i.e., up and down 
the array). Its shape differs slightly in the two nod 
positions because the nod positions are at differ- 
ent distances from the center of the focal plane. 
The extraction algorithm interpolates or extrapo- 
lates the trace polynomial coefficients depending 
on the source position. Figure[2] presents the PSF 
for each module. Our reconstruction of the PSF 
is good enough to reveal multiple Airy rings. Fig- 
ure^ displays cuts along various wavelength rows. 

3. Methodology 
3.1. Wavelength grid 

When performing the spectral extraction, the 
algorithm considers each row independently. The 
wavelength actually varies along a given row so 
that a different location in the spatial direction 
_p <jorresponds to a slightly different wavelength. In 
theory the use of pseudo-rectangles is adequate 
to consider windows with uniform wavelength but 
the practical implementation is limited by the 
lack of precise knowledge of the response function 
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Fig. 3. — Cut of the SL and LL super-sampled PSF along various wavelength rows. 



Table 2 

Nod position coordinates and trace equations. 



Field of view 


Nod position 




Equation 




SL1 nod 1 
SL1 nod 2 


9.26 
19.85 


(4.51 x 1(T 
(4.51 x 10" 


2 )y - (2.21 X 10" 
2 )y - (2.21 X 10" 


5 )y 2 
5 )v 2 


SL2 nod 1 
SL2 nod 2 


55.53 
65.52 


(4.32 x 10" 
(4.32 x 10" 


' 2 )y - (3.71 x 10" 
2 )y - (3.71 x 10" 


7 )v 2 


SL3 nod 1 
SL3 nod 2 


52.34 
62.46 


(4.50 x 10" 
(4.50 x 10" 


2 )y - (1.90 x 10" 
2 )y - (1.90 x 10" 


V 

5 )v 2 


LL1 nod 1 
LL1 nod 2 


36.29 
47.00 


(-5.10 x 10 
(-4.71 x 10 


- 2 )y + (4.25 x 10 
" 2 )y + (3.93 x 10 


-"h 2 

- 5 )v 2 


LL2 nod 1 
LL2 nod 2 


74.02 
85.00 


(-6.69 x 10 
(-6.16 x 10 


- 2 )y - (5.58 x 10 
- 2 )y - (5.14 x 10 


- h )y 2 
- & )v 2 


LL3 nod 1 
LL3 nod 2 


74.09 
85.20 


(-6.69 x 10 
(-6.16 x 10 


- 2 )y - (5.58 x 10 
- 2 )y - (5.14 x 10 


- b h 2 

- 5 )y 2 



NOTE. — Nod positions arc given in pixels and are derived from the trace 
intersect at row 0. y represents the row number in the detector array. 
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within a pixel (see introduction). Furthermore, 
a slight loss of spectral resolution occurs due to 
the artificial splitting of information in individ- 
ual pixels between wavelength elements. Regard- 
less of these side-effects, considering rows instead 
of pseudo-rectangles provides a valuable means to 
better sample the resulting spectrum by using the 
exact wavelength at the nod positions. 

Figure^] shows an example of a spectrum re- 
sulting from the combination of two nod spectra 
in which the wavelength is determined directly at 
the source location in the cross-dispersion direc- 
tion. The slight wavelength shift between the two 
nod positions provides a better spectral sampling 
at the expense of a lower S/N ratio because the 
nod spectra are not combined in the same wave- 
length scale. Narrow spectral lines benefit greatly 
from the full sampling as their profiles are bet- 
ter determined, implying a better detection level, 
and a more accurate line flux measurement. It 
is necessary that the two nod spectra align fairly 
well for them to be combined. For this reason, we 
implemented an automatic algorithm which, pro- 
vided two nod spectra, fits the error function (dif- 
ference between the spectra) with a 6-spline. The 
interpolated difference is then split evenly to the 
individual spectra which can then be combined. 
This algorithm preserves the spectral resolution. 
Practically, the SMART software can regrid a poste- 




Wavelength [\xm] Wavelength \\im] 

Fig. 4. — SL spectrum of the galaxy NGC 6240 re- 
sulting from the combination of the two nod spec- 
tra. The blended lines H2 and [Ar 11] are shown on 
the left while the PAH feature at 11.3 /im is shown 
on the right. 



riori the data on the reference wavelength grid or 
keep instead the wavelength determination at the 
source position (see online documentation). 

3.2. Optimal extraction 

Before the PSF can be used to perform the op- 
timal extraction, it has to be downsampled and 
shifted to the source position. This important 
process consists of three steps: (1) the data grid 
(sub-array of the original 128-pixel row where the 
spectral order is the one requested) is stretched to 
match the (super-sampled) PSF grid and is shifted 
to align the PSF and the source. (2) The PSF is 
interpolated in the new grid. (3) The PSF flux 
of a given pixel in the data grid is determined by 
integrating the flux of the subpixels. 

To find the flux from the source in each row, 
we use a multiple linear regression algorithm. The 
scaling factor /(A) gives the flux for a given row, 
and is thus constant along the row. Each pixel can 
thus be used to determine /(A) independently: 

/(A) = Di/Pi Vi e [0,n], (3) 

where i is the column element of the slit spatial 
profile, D is the source profile, and P is the PSF 
profile. The scaling factor /(A) provides the source 
flux in e~/sec, which eventually leads to the source 
flux density in Jy after calibration. The flux cal- 
ibration performed by AdOpt uses the calibration 
files b{0 2}_aploss_f luxcon . tbl provided by the 
SSC. These files provide a calibration which does 
not include extraction aperture-loss light correc- 
tions due to the size of the regular extraction win- 
dow. It is thus well suited for the purpose of op- 
timal extraction which does not involve a regular 
aperture. An additional relative spectral response 
function (RSRF) is used to correct for the small 
(-10%) residual flux offsets. The RSRF used in 
SMART is calculated by performing optimal extrac- 
tion on the calibrator star HR 7341 which is then 
compared to its stellar theoretical template pro- 
vided by the SS(M (see also Decin et al. 2004). 
The flux calibration may be subject to slight mod- 
ifications in the future as new versions of calibra- 
tion files will be provided by the SSC. 

The use of optimal extraction increases the 
overall quality of the spectra, with the most signif- 
icant improvements for observations with low S/N. 

5 http://ssc.spitzer.caltech.edu/IRS/calib/templ/ 
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Fig. 5. — Illustration of the S/N improvement provided by optimal extraction. The object is the galaxy 
NOA014213. 5+35.1. The histogram shows the tapered column extraction while the solid line spectrum 
shows the optimal extraction. The improvement in the S/N ratio is a factor ~ 1.8 in this case. 



Quantitatively, the S /N ratio increases by a factor 
~ 1.5 — 2.0 for sources as faint as 1 mJy compared 
to normal tapered column extraction. In addition, 
while tapered column extraction is sensitive to bad 
pixels, optimal extraction is mostly unaffected, ex- 
cept when the bad pixels are precisely located at 
the peak of the spatial flux distribution. 

3.3. Multi-source extraction 

One of the most interesting applications en- 
abled by the knowledge of the super-sampled PSF 
is the possibility of extracting spatially-blended 
sources. In the case of multiple sources, the es- 
timated spatial profile is given by 



D, 



k 



A (A) x P s 



(k) 



(4) 



where m is the number of sources, fk is the fit pa- 
rameter, and P( fc ) is the PSF aligned with the k-th 
spectrum. As explained by Collins, Gull, Bow- 
ers & Lindler (2002) and Bcvington & Robinson 
(1992), the coefficients can be found by solv- 
ing the multiple linear regression using the least- 
square method: 

/ = (PP T )- 1 P T D, (5) 

where D is the vector {Di} and P is the ma- 
trix {P^}. The syste m is characterized by n 



equations (number of columns where the order 
is as requested) for m unknowns. The system 
thus remains overdetermined for a large number 
of sources. We find that it is possible to ex- 
tract spectra of sources separated by two pixels 
or more along the spatial axis. Sources separated 
by one pixel can be extracted, with recognizable 
spectral features, but the measurement of the fea- 
tures (integrated flux, equivalent width) should be 
regarded with caution. As a result of the min- 
imum separation between two source and of the 
system overdetermination, the algorithm can ex- 
tract a number of sources which is less than half 
the number of pixels in the row (minus the num- 
ber of degrees of freedom for the background, see 
Sect. 13.5]) . Thus, theoretically the maximum num- 
ber of sources that can be extracted is ~ 17 for SL 
to « 15 for LL. Figure[6] gives an example of a real- 
case scenario in the Small Magellanic Cloud, while 
Figure[7] represents an artificial case for which the 
recovered spectra are similar to the original spec- 
tra used to construct the artificial test image. 

There are generally two main scenarios adapted 
to source deblcnding. 

1 . Sources can be known a priori in which case 
their position on the sky can be used to 
extract their spectra. Sources can be also 
identified in the detector image when spatial 
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Fig. 6. — Real-case scenario of slightly blended 
sources in the Small Magellanic Cloud. The spa- 
tial profile is calculated based on the collapse of 
the shortest wavelengths, i.e., where the PSF is 
the narrowest. 
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Fig. 7. — Spectra extracted from an artificial im- 
age in which spectral images of the quasar 3C371 
and the planetary nebula SMP-LMC35 have been 
combined after shifting one arbitrarily by two pix- 
els in the cross-dispersion direction. The resulting 
spatial profile is shown on the right. 
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structure is evident. Note however that the 
multi-source extraction is originally meant 
only for point-like sources. 

2. Sources can be discovered with the help of 
the source finder (residual minimization, see 
Sect.[3j|. 

Because the multi-source extraction is an in- 
teresting tool for analysis of crowded regions, the 
AdOpt software allows the user to display the slit 
projected on the sky along with the field of view 
position and the requested coordinates. Simbad3 
and NEE0 sources in the field can also be displayed 
for easy source identification. Finally, any FITS 
image can be overplotted, including automatic 
queries from the SDSS, 2MASS, Spitzer/IRAC, 
and Spitzer/MIPS catalogs. 

3.4. Spectral-pointing induced through- 
put error (SPITE) 

No calibration currently exists for offsets in the 
dispersion direction. Such offsets can result in sig- 
nificant light loss in the SL module when sources 
are more than ~ 1" from the slit center (see IRS 
reports by Sloan 2004; Sloan, Nerenberg & Rus- 
sell 2003; Nerenberg & Sloan 2003). This effect is 
known as the spectral-pointing induced through- 
put error (SPITE). Since the LL aperture is much 
wider than the SL one (TableQ]), the effect is unim- 
portant for LL except for significant offsets. 

Accounting for offsets in the dispersion direc- 
tion keeps the optimal extraction valid for sources 
located anywhere in the aperture. It provides a 
reliable flux calibration which is useful in the case 
of a mispointing, and it is also adapted to the case 
of multiple sources as it is likely that at least one 
source is not near the slit center. While a regular 
tapered extraction would require only the knowl- 
edge of the flux fraction that is lost outside the slit 
(provided by the b{0 2}_slitloss_convert .tbl 
tables from the SSC), optimal extraction requires 
the knowledge of the PSF profile at the exact 
source position. Available data could not be used 
to derive the PSF profile at various positions in the 
dispersion direction. We used instead a theoreti- 
cal PSF and calculated the function transforming 
the "nominal PSF" (i.e., integrated PSF across 

6 Simbad is found at http://simbad.u-strasbg.fr/simbad/ 
7 NED is found at http://nedwww.ipac.caltech.edu/ 



the slit for a source at the slit center) to the effec- 
tive PSF (source anywhere across the slit). The 
function is then applied to the empirical super- 
sampled PSF derived (Sect[2]). Finally, the effec- 
tive PSF is normalized to the total area of the 
nominal PSF so that the scaling factor provides 
the correct flux value without any further calibra- 
tion required other than the one already applied. 

With the throughput error corrected, it be- 
comes possible to extract sources anywhere in the 
slit with a reliable flux calibration. SMART allows 
the user to extract sources at specific celestial co- 
ordinates, but can also find the source in the dis- 
persion direction. This is possible because as the 
source moves away from the slit center, the effec- 
tive PSF becomes wider so that the quality of the 
fit varies. Figure[5] shows an example on an actual 
mispointing of sal. 6" in the dispersion direction of 
the calibration star HR6348. The offset is about 
half the SL aperture height, providing an adequate 
illustration of the algorithm explained above. By 
using an on-thc-fly PSF, the flux level agrees with 
the true spectrum while providing smaller residu- 
als (by 20%). 

3.5. Complex background 

There are usually several ways to remove the 
background. The background contribution can be 
inferred by interpolating the columns on each side 
of the source, which requires good knowledge of 
the source position. Another way to remove the 
sky relies on the subtraction of dedicated back- 
ground observations or differencing the image by 
nod or by order (e.g. SL nod 1 minus SL nod 2, 
or SL1 nod 1 minus SL2 nod 1). The disadvan- 
tage of these methods is that they assume that 
the background at the position where it is mea- 
sured is the same to where the source is located. 
The new algorithm can calculate the underlying 
extended emission together with the PSF scaling 
factor (see also Geers et al. 2006, Lahuis 2007). 
Therefore, even the pixels where the source emits 
are taken into account for constraining the back- 
ground shape and level. 

The AdOpt algorithm assumes that every col- 
umn in a given spatial row should follow the same 
background equation. Any background shape can 
be considered which is a function of the column 
index i, i.e., with £>j = a/ x g(i, I) with the co- 
efficients a; being the unknowns and g(i, I) repre- 
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sents any function that describes the background. 
We thus have again n linear equations: 



A = X>(A) xPf'+^a, xg(i,l). (6) 




Wavelength [um] 




140 160 1 3D 200 

Column index (subpixel space) 



Fig. 8.— The SL spectrum of the star HR6348 is 
shown on the left. The reference spectrum shows 
the median spectrum of all the existing observa- 
tions of this star. The "Centered PSF" spectrum 
shows the extraction of the mispointed exposure 
assuming that the PSF is centered in the disper- 
sion direction. The "Effective PSF" shows the ex- 
traction of the same image, by calculating the PSF 
on-the-fly at the real source position. The corre- 
sponding PSF profiles are shown on the right. 



The system is solved using the multi-linear regres- 
sion algorithm, which estimates simultaneously 
the unknowns fk and a/. For practical reasons, 
the algorithm AdOpt makes use of a polynomial 
background, i.e., g(i,l) = i . The system remains 
overdetermined if the sum of the number of sources 
and the background parameters does not exceed 
w 17 for SL to w 15 for LL. An example of an 
optimal extraction of a source within a complex 
background is shown in Fig. [9] 

The particularity of this method is that the 
background profile is estimated for each row in- 
dependently, which can lead to undesirable small- 
scale variations along the wavelength axis. Hence 
for noisy images, it is advised to perform a first ex- 
traction to estimate the row-by-row background 
and then smooth the reconstructed background 
image along the wavelength axis (see online docu- 
mentation) . 

3.6. Source finder 

The plug-in AdOpt includes a new source- 
finding algorithm. Several source finders already 
exist, using for instance the centroid in the col- 
lapsed spatial profile (along the PSF trace), or the 
median of the centroids in each row. These meth- 
ods require fitting a profile, usually using a Gaus- 
sian curve, which inevitably introduces possible 
systematic errors and even biases. These biases 
can introduce negative effects in the data, such 
as noticeable bumps in the spectrum (which are 
more abrupt the larger the discrepancy between 
the source finder and the source actual position) 
(see example in Figure fTU)) . 

The new source finder uses the super-sampled 
PSF. The PSF is aligned at varying positions along 
the slit. For each position, the residual image 
(data image minus the PSF image) is then ana- 
lyzed. First, the residual image is rectified using 
the trace equations of Tabled Then, a &-spline is 
fitted to each column, with rejection parameters 
such that bad pixels are identified and ignored. 
The median value of the fit is then calculated for 
each column, which results in a representative spa- 



10 



coordinate wrt roc position 




-coordinate (@naw 0) 



Fig. 9. — Real-case example of how a complex 
background is handled by SMART. The background 
is determined together with the optimal extraction 
of the sources (two stars in the Small Magellanic 
Cloud). 
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Fig. 10. — Spectrum of the quasar 3CQ-3C371 re- 
sulting from the extraction of a single DCE image. 
Top — The dip at short wavelengths observed is 
typical of the artifacts introduced by slight uncer- 
tainties in the source position coupled with a nor- 
mally sampled PSF. Bottom — The artifact could 
be exactly reproduced by using the super-sampled 
PSF of SMART and a forced shift of 0.1 pixel. 



tial row. This row is ideally zero when there is no 
background and when the PSF has been aligned 
to the source. Hence, by simply taking the aver- 
age of this row, one has a very good leverage on 
the quality of the alignment PSF-sourcc. The pre- 
cision reached by the source finder is better than 
l/20 th of a pixel. 

4. Other applications 

Besides the improved S /N and the multi-source 
extraction, the precise knowledge of the PSF pro- 
file enables several valuable possibilities for the 
handling of point-like sources while providing in- 
teresting tools for handling extended sources. 

4.1. Spatial extent 

With the knowledge of the super-sampled PSF, 
the spatial extent of a given source can be de- 
termined fairly accurately by comparing its spa- 
tial profile to the PSF for the same wavelength or 
wavelength range. The plugin AdOpt allows the 
user to check the extent of any source within the 
slit, by providing three meaningful quantities: the 
percentage of the extent with respect to the PSF, 
the extent in pixels, and the extent in arcseconds. 
Among the most useful applications to check the 
source extent is to identify the possible presence of 
another source or the presence of extended emis- 
sion associated with the point-like source (sec 14. 21) . 

4.2. Point-like sources embedded within 
an extended source 

The complex case of point-like sources em- 
bedded within an extended emission can be ap- 
proached in various ways depending on the back- 
ground nature. 

• Large-scale background with spectral fea- 
tures. A low-grade polynomial background 
can be estimated in every wavelength row 
(Sect. 13.5) ). The image quality should how- 
ever be better than some threshold. 

• Large-scale background with no spectral fea- 
tures. The same procedure as above is ap- 
plied, but the final background image is 
calculated by smoothing the "row-by-row" 
background. 
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• Small-scale background (e.g., disk, extended 
region, ...). A polynomial of high degree 
should be used, implying that less con- 
straints will be possible. In fact, a polyno- 
mial might not suffice to adequatly represent 
the background, but it will still allow a first 
order approximation (see next section). 

4.3. Removal of point-like source contri- 
bution 

In the case of a point-like source embedded in 
an extended source, it is possible to extract the 
spectrum of the extended source alone by extract- 
ing first the point-like source (no background sub- 
traction in this case). The reconstructed image 
(which in the case of images with normal point 
sources contains the residuals) includes only the 
extended source. Note that the spectrum of the 
point-like source includes a fraction of the ex- 
tended emission and should be regarded with cau- 
tion, especially for flux calibration and spectral 
features measurements. The extended emission in 
the reconstructed image can then be extracted us- 
ing either a full slit or a tapered column. It is 
important to bear in mind that there is no flux 
calibration available for the extended emission, ex- 
cept if the source fills uniformly the aperture. 

5. Summary 

We have presented the optimal extraction al- 
gorithm adapted to the Infrared Spectrograph on- 
board Spitzer. Besides providing a significant in- 
crease in the signal-to-noise ratio as compared to 
the regular "tapered" extraction, the new algo- 
rithm allows extraction of multiple sources at any 
location within the aperture. The optimal extrac- 
tion is also implemented for sources shifted in the 
dispersion direction. Finally, it is possible to per- 
form an optimal extraction with a complex ex- 
tended background emission. 

The AdOpt plugin is released with the SMART 
package (versions equal or later than 8.0). The 
code is available at the Infrared Science Cen- 
ter website, along with extensive documentation. 
SMART and the optimal extraction will be main- 
tained in the future, with a possible inclusion 
of the optimal extraction for the high-resolution 
modules. 
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